Singlet Generation in Mixed State Quantum Networks 
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We study the generation of singlets in quantum networks with nodes initially sharing a finite 
number of partially entangled bipartite mixed states. We prove that singlets between arbitrary 
nodes in such networks can be created if and only if the initial states connecting the nodes have a 
particular form. We then generalize the method of entanglement percolation, previously developed 
for pure states, to mixed states of this form. As part of this, we find and compare different distillation 
protocols necessary to convert groups of mixed states shared between neighboring nodes of the 
network into singlets. In addition, we discuss protocols that only rely on local rules for the efficient 
connection of two remote nodes in the network via entanglement swapping. Further improvements 
of the success probability of singlet generation are developed by using particular forms of 'quantum 
preprocessing' on the network. This includes generalized forms of entanglement swapping and we 
show how such strategies can be embedded in regular and hierarchical quantum networks. 

PACS numbers: 03.67.Bg, 03.67.-a, 64.60.ah 



I. INTRODUCTION 

Quantum entanglement is one of the most notable fea- 
tures of quantum systems and has been accepted as a key 
resource for quantum information processing [l[ . The dis- 
tribution of entanglement through quantum networks is 
therefore essential for the future of a variety of applica- 
tions, ranging from quantum cryptography to quantum 
teleportation and distributed quantum computing 
However, the generation of these entangled states faces 
a severe obstacle. Quantum channels such as free-space 
transmission or optical fibers are prone to loss and de- 
coherence. This causes the desired maximally entangled 
states to degrade into mixtures and limits the distance 
over which the quantum information can be sent directly. 
To overcome these problems 'quantum repeater' schemes 
have been proposed [1, 3, M , 0, S Hi which make use of 
the ability to 'purify' |9l.]lol| and 'swap' [O, entangle- 
ment to maintain a high fidelity throughout. Quantum 
repeaters are a promising tool for entanglement distri- 
bution, particularly since the amount of required phys- 
ical resources increases only polynomially with the dis- 
tance @, but operate in a ID setup of network nodes. 
Real networks are typically two- (or higher) dimensional 
and it is therefore desirable to study if entanglement dis- 
tribution can be made more efficient in these cases. 

A scheme for entanglement distribution in higher di- 
mensional networks was recently proposed by Acm et 
al. [l^ in which ideas from classical bond percolation 
have been applied to regular, i.e. lattice-shaped, quan- 
tum networks. The scheme makes use of the networks' 
connectivity and allows for the generation of maximally 
entangled singlet states between arbitrary points of the 
network, with a probability that is independent of their 
separation. The only requirement is that the nodes are 
initially connected by bipartite pure states with suffi- 
ciently high entanglement. The restriction to pure states 
was made since a pure, partially entangled state can be 



converted into a singlet with finite probability via lo- 
cal operations and classical communication (LOCC) [l3| 
which is essential for the bond percolation protocol: Ini- 
tially one attempts to convert all bipartite pure states 
into a singlet which, in each case, succeeds with a cer- 
tain probability. If this singlet conversion probability 
(SCP) exceeds a lattice-geometry-dependent threshold, 
arbitrarily large clusters of singlet-connected nodes form 
which can successively be connected via entanglement 
swapping. In this way we can create a singlet between 
arbitrarily remote nodes in the network. However, it was 
pointed out in [Tsj that this process, known as Classical 
Entanglement Percolation (CEP), is not optimal since 
certain quantum preprocessing schemes a pplie d to the 
network can improve the SCP [H, [H, [H, [11, 

and thus it is possible to apply bond percolation to lat- 
tices in which this would otherwise not be possible. 

Clearly, the assumption of having a pure-state net- 
work is an idealization and in any practical situation the 
states connecting the nodes of the network will be mixed. 
In the idea of entanglement percolation was applied 
to mixed states for the first time. In this paper we elab- 
orate and extend the ideas presented in [20j. The net- 
works we consider are composed of nodes, each of which 
can consist of several qubits, and may be connected by 
a finite number of bipartite mixed states (see Fig. [1]). 
We aim to create a perfect singlet between two arbitrary 
nodes in the network using a finite amount of resources, 
i.e. a finite number of initial states which are converted 
into a singlet, which distinguishes our and other entan- 
glement percolation schemes from, e.g., the quantum re- 
peater protocol where one aims to generate a state with 
high but non-unit entanglement fidelity. Particularly we 
structure the paper as follows. 

In Sec. [H] we prove a necessary and sufficient condition 
that a perfect singlet can be generated in a network of 
arbitrary geometry the nodes of which are initially con- 
nected by bipartite mixed qubit states. We show that 
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singlet generation between two nodes is possible if and 
only i/they are connected by at least two 'paths' consist- 
ing of a particular class of states. These states arise nat- 
urally in systems undergoing amplitude damping. Thus 
our result is not only of theoretical but also of practical 
relevance. Unfortunately, the proof does not deliver an 
efficient scheme for singlet generation. We therefore spe- 
cialize in the remaining sections on networks with regular 
geometry, i.e. lattices in 2D and 3D and devise general- 
izations of entanglement percolation to the mixed states 
described in Sec. [ITl 

In Sec. mil we briefly summarize the idea of classical 
entanglement percolation with pure states. 

In Sec. llVl we extend the concept of classical entangle- 
ment percolation to mixed states. To this end we con- 
sider regular networks where each node is connected to its 
neighboring nodes by a finite number of the mixed states 
introduced in Sec. [TTl We present two different distilla- 
tion protocols which are used to convert these states into 
a singlet with a probability above the percolation thresh- 
old of a variety of lattice geometries. After the distilla- 
tion, clusters of singlet-connected nodes emerge and we 
aim to create a singlet between two nodes in such a clus- 
ter by successive application of entanglement swapping. 
By communicating classically each node can determine if 
singlets exist between it and its neighboring nodes. This 
information can be communicated and stored classically 
in a central data processor. Typically one would then use 
this information to apply a path-finding algorithm which 
locates a suitable 'path' of singlets before swapping oper- 
ations are performed. As an alternative to this we discuss 
a classical and a quantum protocol which merely require 
classical communication between neighboring nodes and 
basic computing within each node. The quantum proto- 
col relies on the formation of many-qubit GHZ states via 
local operations and classical communication with neigh- 
boring nodes and subsequent measurements at all nodes 
except the ones to be left in the final singlet. 

In Sec. |V] we show that the idea of 'quantum pre- 
processing' as it was successfully applied in pure state 
networks can be generalized to mixed states. In par- 
ticular we devise a number of strategies on small net- 
works which improve the SCP, and we show that these 
smaller networks can be embedded into larger networks 
to enable CEP which would otherwise not be possible. 
Furthermore, we discuss 'hierarchical schemes', i.e. net- 
works which are defined iteratively and were first dis- 
cussed in [l^, [l^ . Also in these cases it turns out that 
quantum methods outperform classical percolation. Fi- 
nally, in Sec. IVII we summarize and conclude. 



II. SINGLET GENERATION WITHIN AN 
ARBITRARY MIXED STATE NETWORK 

In this section we consider quantum networks of arbi- 
trary geometry as shown in Fig. [T] where the qubits in the 
nodes are 'connected' by bipartite mixed states to qubits 




FIG. 1: Mixed-state quantum network. Qubits in a node 
(circles) may be connected by bonds (thick lines), i.e. they 
share mixed entangled states, 'edges', (solid, black lines) of 
qubits (black dots) with other nodes. When two 'paths' of 
states of the form ([l} connect A and B a singlet (dashed 
line) can be created with finite probability. This is proven by 
partitioning the nodes into two groups with one containing 
A (shaded region) and the other B. For it to be possible to 
generate a singlet between A and B these groups must be 
linked by at least two states of the form ^ for all possible 
partitions. 

in other nodes. We will call a single bipartite mixed state 
an edge and the set of edges directly connecting two nodes 
a bond. Note that an edge connects exactly two qubits 
in different nodes. In the following we will prove that 
the generation of a perfect singlet between two arbitrary 
nodes A and B with finite probability in such a network 
is possible if and only if there are at least two paths of 
states linking A and B which have, up to local unitaries, 
the form 

p(a,7,A) = A|a,7)(a,7| + (1 - A)|01)(01|, (1) 

where \a,-f) = y^\OQ) + y/l - a- 7|11)+^|01) and < 
A < 1. We show this by separately proving a necessary 
and sufficient condition which, together, prove the above 
statement. 

Necessary condition. We split the network into two 
groups of nodes. A, containing A and a finite number 
of other nodes, and B, which consists of the rest of the 
network and particularly contains B. These groups are 
linked by a finite number of edges. A singlet can be 
established with finite probability, via local operations in 
the groups and classical communication between them, if 
and only if at least two of the states have the form ([T|). 
Appendix |X] contains a concise proof of this fact based 
on [Hj which agrees with the result of Ref. [l^, that. 



in general, a singlet can not be generated with a finite 
probability from a finite number of mixed states. 

With two states of the form p{a,j,X) and 

p{P,S,h'), we obtain a singlet with a finite probability 
by first performing two C-NOT gates locally, with the 
p{/3, 6, ly) state's qubits acting as the target qubits. These 
target qubits arc then measured in the computational ba- 
sis. If wc find both qubits to be in the state |1) we have 
generated a pure entangled state between the qubits that 
originally corresponded to the p{a, 7, A) state. We will 
refer to this measurement as the pure state conversion 
measurement (PCM). The state formed is 



\a') = \a'r/ = 0) = V^\oo) +VT~^\n), (2) 

i.e. a' is a Schmidt-coefficient that has the value 
, _ min(a(l - ^ - (5), /3(1 - a - 7)) 



a(l - /? - (5) + /3(1 - a - 7) 



(3) 



The probability that the PCM succeeds in generating this 
state is given by 



Pc 



Ai^(a(l - (3 - S) + f3{l - a - 7)). 



(4) 



For identical states, i.e. a = f3, j = 6, the PCM already 
yields a singlet. Otherwise the state can be transformed 
into a singlet via the 'Procrustean method' [2^ that con- 
verts any pure 2-qubit state \a') into a singlet |l/2) with 
a probability 2min(l — a', a'). The total success proba- 
bility of generating a singlet is then given by the SCP 



Pc 



2Ai^min[a(l-/3-5),/?(l-a-7)] (5) 



which coincides with the optimal probability for creating 
a singlet from two of these states (2]| . 

We can perform this partition of the network in an 
arbitrary way, as long as one group contains A and the 
other contains B. To be able to create a singlet between 
A and B via LOCC we must have at least two states of 
the form ([T]) in all possible partitions. This gives us a 
necessary condition that to create a singlet between two 
nodes with a non-zero probability there have to be at 
least two distinct 'paths' of edges of the form ([T|) con- 
necting the corresponding nodes. In Fig. [Ija) this is in- 
dicated by two spatially distinct paths of bonds. The 
states of the qubits that are not contained in this path 
are irrelevant and can therefore be in arbitrary states. 

Sufficient condition. In order to show this we make 
use of entanglement swapping. This operation can be 
performed in the setup shown in Fig. [2] and consists of 
performing a measurement in the standard Bell basis on 
the qubits located at C2 and LOCC which causes Ci 
and C3 to become entangled. If the edges are of the 
form (IT|), p(q;,7. A) and p{(3,5,v), then there are four 
possible outcomes. The probabilities to obtain measure- 
ment outcomes corresponding to the Bell states l^*^) = 
(|00) ± |ll))/%/2 and |$±) = (|01) ± |10))/\/2 are 

P(*±) = \{h±\v+{l-P-5){l-\)y + aX{l^v)) (6) 



Ci C2 C3 



FIG. 2: Basic arrangement for entanglement swapping. En- 
tanglement swapping involves a measurement in the Bell basis 
at node C2 and classical communication between the nodes 
followed by local unitaries which causes Ci and C3 to become 
entangled. 



and 



P($±) = -{9±\y + {l-v)[l-aX) + {l3 + 5){l-\)u), (7) 



where 



h± =a(i + (1 - a - 7)(1 ~P-5) 



(8) 



9± ^Jf3 + (1 - a - 7)'5 + (1 - a - 7)/3 

+ i^±^ail-f3~6))-\ (9) 

If we measure the qubits at B to be in the states l^**) 
then we actually form another state. 



^(3 (V^±^7(l-/3-(5))2 Xiyh± 



(10) 



of the form ^ between Ci and C3. Unfortunately for 
the other outcomes the states' form is not generally main- 
tained. Note that if S ~ j — we can discard these cases 
by replacing the state with |01) leading to an operation 
that transforms p{a, 0, A) (g) p{P, 0, v) into 



'a/3 ^ , 
P I j^,0,Xiyh± 



(11) 



which will be useful in Sec. IV Al We can therefore create 
a state of the form ([T]) with non-zero probability between 
two nodes of the network, e.g. A and B in Fig. [T] given 
that these nodes are connected by a path consisting of 
states of the same form. Two such states, originating 
from two paths, can then be converted into a singlet, 
using a PCM and the Procrustean procedure. Unfor- 
tunately, this scheme leads to an exponential decrease 
of entanglement fidelity and thus success probabil- 
ity, with the number of swapping operations. Hence it is 
not an effective solution to the problem of long-distance 
entanglement distribution. In Sec. IIVI we will therefore 
introduce effective protocols which can be applied in reg- 
ular network geometries and succeed in creating a singlet 
with a probability independent of distance. 

Note that when entanglement swapping is done with 
pure states all of the outcomes can be used, and if these 
outcomes occur with probabilities Pm the pure state \a) 
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FIG. 3: Illustration of classical entanglement percolation in 
a square network. Pairs of qubits (black dots) in neighbor- 
ing nodes (circles) are in identical, pure, partially entangled 
states (solid, black lines). The percolation scheme involves 
these entangled states being converted into singlets (dashed 
lines) with probability p. If p exceeds the percolation thresh- 
old these form large clusters and we can obtain a singlet be- 
tween any two qubits within a cluster by performing swapping 
operations. 



with 



is recovered by using classical communication and local 
unitaries. 



any two nodes in the cluster. In the theoretical case 
of an infinitely large lattice a cluster that is infinite in 
extent forms if and only if p > pc, where Pc is a lattice- 
dependent percolation threshold. This approximates the 
case for large but finite lattices where the threshold be- 
comes more definitive as the size of the lattice increases. 
Values of pc for a number of lattice geometries are given 
in Table HI If each bond in a network consists of a single 
pure state ja) we can calculate a threshold for a given 
by 2 min(Q;, 1 — a) > Pc- The probability that a node be- 
longs to the infinite cluster is known as the percolation 
probability Two randomly chosen nodes are both 

part of the infinite cluster with a probability 0(p)^ and 
thus can be connected over an arbitrary distance. 



Lattice 



Threshold pc 



2D Square 
2D Triangular 
2D Honeycomb 
3D Simple Cubic 



0.5 

2sin(7r/18) ^ 0.347 
1 - 2sin(7r/18) 0.653 
0.249 



3D Face-Centered Cubic ^ 0.120 

TABLE I: Threshold probabilities for various regular network 
geometries [23.[25t. 

It has been shown that CEP using pure states is not 
optimal and that by performing particular quantum pre- 
processing steps, particularly swapping operations on 
the lattice before converting to singlets, improvements 
can be achieved. These improvements include obtain- 
ing a geometry with a lower percolation threshold af- 
ter the swapping operation and splitting the lattice into 
two, so that a higher percolation probability can be ob- 
tained [H, [l^, [Ty, [13, [l^l- Recently, another method, 
that transforms the initial bipartite network into a prob- 
abilistic multipartite network, has also been shown to 
yield an improvement p^ . 



III. CLASSICAL ENTANGLEMENT 
PERCOLATION WITH PURE STATES 



IV. CLASSICAL ENTANGLEMENT 
PERCOLATION WITH MIXED STATES 



In this section we will briefly review the use of perco- 
lation for distributing singlets in pure state networks [l^, 
[isj . known as classical entanglement percolation (CEP). 
The procedure is based on classical bond percolation, 
where we consider a regular lattice of nodes connected 
by identical quantum states, as shown in Fig. [3l A de- 
scription of classical bond percolation can be found in 
Ref [24]. If the nodes arc connected by pure states of 
the form |a) they can be converted into singlets using 
the Procrustean method with a SCP p = 2 min(a, 1 — a). 
These sing lets act as the bond in the bond percolation 
model |38| and are distributed randomly with a proba- 
bility p. The nodes that can be connected by a path 
of singlets form a cluster. By using entanglement swap- 
ping (sec Sec. |ll| we can then generate a singlet between 



In this section we extend CEP to mixed states. We 
consider regular lattices, e.g. triangular (see Fig. d]), 
square, or even lattices in higher dimensions. Bonds 
between network nodes are composed of multiple edges 
which satisfies the necessary condition proven in Sec. [TTl 
We assume that each bond is identical. When these 
bonds contain at least two states of the form ([T|) they 
can be converted into singlets by PCM followed by the 
Procrustean method. If the probability that a bond be- 
comes a singlet exceeds the percolation threshold CEP is 
achieved. In the remainder of the paper we will assume 
that the states forming edges are of the form ([1]) with 
7 = 0. Setting 7 = is not a major restriction but al- 
lows us to keep the equations manageable. All protocols 
presented in this paper can also be performed if 7 7^ 0. 



5 




1.0 



O 0.6 
m 



0.2 



; Honeycomb 




Sc[iiare 
Triangular 



FIG. 4: Triangular network. This is a simple 2D arrangement 
of PMSs in which CEP is possible. 



FIG. 5; The purification setup consists of n PMSs (solid lines) 
shared between two nodes A and B. The aim is to distill these 
states into a singlet. 



We will call states of the form ([T]) with 7 = 0, i.e. 

p(«,A) = p(a,7 = 0,A), (13) 

purifiable mixed states (PMSs). Note that these states 
form the states of two entangled atomic ensembles in the 
DLCZ quantum repeater scheme 



FIG. 6: Singlet conversion probability for n-edged bonds 
using the recycling scheme for ol = 1/2 and n = 
2, (3),4,6,8, 10, 12, 14, 16 (bottom to top). The n = 3 line 
(dashed) corresponds to the DSS scheme. The percolation 
thresholds for triangular (T), square (S) and honeycomb (H) 
lattices are given by the horizontal lines. 



are 



A^: a|0)l|0)^ + V«(1-«)|1>a|1)b + 



A(l-A) 
A(l-A) 



v/^(r^|2)l|2)^ + (l-a)|3)l|3)|, 
^|0)l|2)| + ^/T^|l)l|3)i, 
V^|0)l|l)i + VT^|2)l|3)i, 

WaWb- (14) 



If this is acted on by the projective measurement 



+ |2)1 



21 at A and 



(1| + |2)|(2| at B the 
state remaining is (|1)1|1)b + |2)^|2)|)/v^. Both of 
these projective measurements only occur with probabil- 
ity 



Vn 



(15) 



Distillation Procedures 



Distillable Subspace Scheme 



We assume that each pair of neighboring nodes is con- 
nected by n PMSs and our aim is to distill these into 
a singlet. The basic setup is shown in Fig. [S] To ac- 
complish this we will use ideas proposed in Rcf. [26l |. 
Here the concept of a distillable subspace (DSS) is in- 
troduced as a subspace such that the local projection of 
the system state into this space is pure and entangled. 
Locating the DSS involves calculating the eigenvectors 
of the state with non-zero eigenvalues. To simplify no- 
tation we will represent the states at A and B using 
the decimal value of its binary form, i.e. for example 
|00110)a|01001)b = |6)1|9)^. 

As an example, in the case of n = 2 identical states 
p{a, A) the eigenvalues and corresponding eigenvectors 



Note that this is the same SCP as obtained for PCM [see 
Eq. ([3])]. For this example there is no choice between 
entangled states to project out and if the original states 
are the same a maximally entangled state is automati- 
cally obtained. For states that arc not identical this does 
not need to be the case. 

An extension of this scheme to n identical copies of 
PMSs p(a. A) yields the SCP 



Pn = Y.x'^-^ii-xy 

m=0 



^n-m-/c(^ _ (" fe") ( ("7") 



E 

fe=i 



-1) 



(;:) - 1 



(16) 



A derivation of this formula is given in Appendix [B) As a 
particular example it is worthwhile to discuss the case of 
three states in more detail. In this case the measurement 
at A is given by a Positive Operator Valued Measure 
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FIG. 7: The recycling scheme consists of splitting the states 
into pairs that are then purified. If no singlets are successfully 
produced some of the states may still have been transformed 
into PMSs and given enough of these the process can be re- 
peated. 



(POVM) with the elements 



(|0>1(0|- 


h|7)l(7 


), 


(ll>l(ll- 


H2>1(2 


)/2, 


(ll>l(ll" 


h|4)l(4 


)/2, 


(|2)1(2|- 


h|4)l(4 


)/2, 


(|3)1(3|- 


H5)l(5 


)/2, 


(|3)1(3|- 


H6)l(6 


)/2, 


(|5)1(5|- 


h|6)l(6 


)/2. 



(17) 



The measurement at B then depends on this outcome 
and creates a maximally entangled state with a certain 
probability. The SCP is obtained by setting n = 3 in 
Eq. (jl6p and is given by 



Pn 



■i\^a{l - a). 



(18) 



Comparing this with the n ~ 2 case [cf. Eq. (fl3|) ] shows 
an increase in the success probability which can be seen 
in Fig. [6l where the dashed Une represents the SCP for 
three identical states. 



2. Recycling scheme 

The SCP using the DSS scheme does generally increase 
with increasing n. However, the scheme does not make 
use of the available resources in the best way. Indeed, 
the SCP Pn can be significantly improved by grouping 
n identical PMSs into sets of m and converting each of 
these sets into a singlet. For example for to = 2 we apply 
the PCM as described in Sec.|TT]on pairs of states which 
converts them into singlets with a probability given by 
Eq. (O. If this fails for a given pair we may still find both 
measured qubits in the state |0) and have generated an- 
other PMS. This PMS can then be used again in another 
purification attempt. To be more precise, starting with 
n copies of a state p{a^ A) (with a > 1/2) we apply a 2- 
state purification protocol on groups of two. If no singlet 
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FIG. 8: Success probabilities for the recycling schemes that 
split the states into pairs (dashed lines) and sets of three 
(solid lines). Shown are the SCPs for 6, 9 and 12 initial states 
(bottom to top) for a = 1/2. 



is obtained the procedure is repeated on the remaining 
PMSs as illustrated in Fig. [T] The coefficients for the 
PMSs after k repetitions, when no singlet is created, are 
given by 



Oik 



1 ' 2ak-i 



ALi(l 



2aLi ' 
2ak- 



2«Li) 



1 - 2Afe_i + 2A2_^(1 - ak-i + a^i) ' 



(19) 



(20) 



where — a and Aq = A. For states of the form 
p(aA;, Afc) ® p{ak, Afc) the probability of obtaining a PMS 
is cfc = l-2Xk + 2{l- ak+al)Xl. If the PCM yields two 
qubits that are measured in different states the purifica- 
tion step between the two PMSs has completely failed. 
The probability of this is given by fk = 2Afe(l — A^). The 
probability of not generating a singlet using this recycling 
protocol on n states of the form p{ai, Xi) is then found to 
be 



E 

k=0 



LfJ 
k 



1) 



(21) 



where Fo{i) ~ 1. Consequently, the probability of suc- 
cessfully generating a singlet by applying the procedure 
to n states of the form p{a, A) is 1 — i^„(0) which is cal- 
culated iteratively. Examples are shown in Figs. [B] for 
a = 1/2. 

Obviously, the states do not necessarily need to be split 
into pairs. For example we can separate all of the states 
into sets of three and apply the three-state DSS distil- 
lation. In case of failure this can yield a PMS state as 
well, which can then be used in later distillation steps. 
There are a variety of ways to combine the three-state 
distillation with the two-state recycling scheme. Here we 
concentrate on the straightforward approach which only 
uses the three-state distillation on every level of the recy- 
cling scheme. The results are shown in Fig. [51 As can be 
seen in most cases the two-state recycling scheme has a 
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higher chance of success and because of this we will focus 
on the pairing arrangement in this paper. 

B. Percolation Thresholds 

Using the purification procedures described above we 
can apply CEP, as described in Sec. HH for lattice net- 
works with multi-edged bonds. In most cases it is advan- 
tageous to use the two-state recycling scheme, except for 
n = 3 where the DSS scheme should be used. From Fig. [6] 
it can be seen that the SCP increases with the number 
of edges per bond and this allows for a larger range of 
values for A and a such that CEP is successful. For dou- 
ble edged bonds the optimal probability of generating a 
singlet is given by 

< 2A^a(l - a) < 1/2. (22) 

When the bonds are composed of three edges, i.e. three 
PMSs between nodes, we have 

< SX^ail - a) < 3/4. (23) 

By comparing these ranges to the percolation thresholds 
we see that a basic successful setup is a double bonded 
triangular lattice (see Fig. |4]) . The double bonds can be 
converted to singlets and if the chance of this is larger 
than the percolation threshold an infinite cluster will 
form. A singlet can then be created between any two 
nodes within the cluster. Thus percolation occurs if 

2X^a{l - a) > 2 sin(7r/18) w 0.347. (24) 

However, the singlet conversion probability never exceeds 
1/2 for two states. Therefore we would require more 
states in other geometries. For example, if we have three- 
edged bonds between each neighboring node we can apply 
CEP to a square lattice. This is because there are param- 
eters such that 3A^Q!(1 — a) > 1/2. Analogously, CEP is 
also possible in honeycomb lattices with three edges per 
bond. 



C. Local Processing Strategies 

The process of creating singlets, randomly replacing 
the initial network bonds, can be run if each node can 
only communicate classically with their neighbors. Each 
node then knows if a qubit that it contains is part of a 
singlet after this procedure has finished. This informa- 
tion can be stored classically within a node but after the 
bonds are distilled we are faced by the problem of find- 
ing a set of singlets that connect our requested nodes, 
A and B. 

If all of the singlet generation data is collected by a 
'controller' then an efficient path finding algorithm can 
be applied to determine a suitable 'path' of singlets link- 
ing the nodes. An example of a suitable algorithm would 




FIG. 9: Procedure to join a node's singlets onto a GHZ state. 
Here the GHZ state is represented by a dotted box. Addi- 
tional qubits (black dots) that are a part of the GHZ state 
are linked to the dotted box by a dotted line. A CNOT gate 
and measurement (both are represented by a shaded oval) are 
performed between the qubit already in the GHZ state and 
those that are part of a singlet (dashed lines). Each mea- 
surement outcome needs to be sent (gray arrow) to the other 
singlet qubit to perform a local unitary (shaded square) . This 
extends the GHZ state to include qubits connected by edges 
to the node being attached. Once this has occurred for each 
qubit in a singlet the process is repeated by sending out a 
signal to repeat the step at each node that was linked by a 
singlet. 

be a Dijkstra scheme [l^ such as the A* path finding 
algorithm The path information can then be used 
to instruct the correct nodes to perform swapping. The 
swapping operations are performed in order from node 
A to B, so that the measurement outcomes only need 
to be communicated along the chain, between neighbor- 
ing nodes. However this procedure requires one classical 
computer to have complete knowledge of the network. 
Instead; it is interesting to note that this does not need 
to be the case as there are algorithms which do not re- 
quire any more classical communication than this 'con- 
troller' method, indeed they do not require a central 'con- 
troller' at all. This can be done not only classically but 
also via a quantum algorithm using multipartite entan- 
glement which we will introduce below. 

A classical path-finding method would use a type 
of breadth-first search algorithm called a burning algo- 
rithm [2^. Node A sends a 'burning' signal to its neigh- 
boring nodes connected by singlets. These nodes keep a 
record of where they received the signal from and send 
out an identical signal to the other nodes that they are 
connected to. We say that the node has 'burned'. If it 
has already received a signal from a different node then 
the additional signal is ignored. This continues outwards 
from A, 'burning' the nodes. Once node B receives the 
signal it replies to the node it came from with a 'swap- 
ping' message. This node can then perform a swapping 
operation and send another 'swapping' signal, together 
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with the Bell-measurement outcome, back to the node it 
received a 'burning' signal from. The path can then be 
traced back along the nodes with swapping performed at 
each step until node A is reached. Both A and B can 
determine if the protocol has been successful. However, 
A and B may not be in the same cluster and they do 
not know if the protocol has failed when the network is 
of infinite size. This is not a problem for finite networks, 
containing N nodes, as A and B can time the steps taken 
and if these exceed 2(A^ — 1) they both know they are not 
in the same cluster. 

Note that no extra information actually needs to be 
transmitted. We can combine the burning algorithm with 
the process of transmitting the distillation protocol in- 
formation. For example, in a double edged network of 
identical edges, A can perform her PCM and if |1) is the 
outcome she assumes she has a singlet and sends a burn- 
ing signal to the node that would contain the singlet's 
other qubit. If a node receives this signal it can perform 
its PCM and determine if there is a singlet there. When 
there is and if it is the first instance for the node it should 
record that entry qubit and repeat the process, perform- 
ing a PCM on the remaining qubits and sending signals 
to those with the |1) outcome. Once B receives a signal 
it can check that a singlet has been created with a PCM 
and then send a swapping signal back as before. During 
the swapping, a node can use the Bell-measurement infor- 
mation received to indicate that a swapping is required 
so no explicit 'swapping' signal is required either. All 
of this information transfer would have been necessary 
as well if a controller algorithm would have been used. 
Hence the generation of the singlet can be accomplished 
by defining rules for each node and allowing them to run 
with nearest neighbor classical communication. This is 
fundamentally different to the controller process and has 
made use of parallel computation to find a path that no 
single node has full knowledge of. 

We will now consider an alternative, quantum algo- 
rithm that is based on the burning algorithm and makes 
use of multipartite entanglement in the network. The 
protocol starts after we attempted to convert all bonds 
into singlets and every node has knowledge about its sin- 
glet connections to nearest neighbors. We build up a 
progressively larger multi-qubit GHZ state, defined by 
\GHZ„) = (|0)i...|0)„ + |l)i...|l)„)/V2, spread between 
the 'burned' nodes by adding qubits in each burning step. 
Building up such a state requires joining two GHZ states, 
\GHZn) and \GHZm), to create \GHZn+m-i) (note that 
a singlet equals IGHZ-z))- This is done by performing 
a CNOT gate between a qubit in \GHZn) and a tar- 
get qubit in \GHZm), measuring the target qubit in the 
Z-basis, communicating the measurement result to the 
other qubits in \GHZm) and performing a unitary op- 
eration on them depending on the outcome. Now we 
perform the same process as for the 'burning algorithm', 
however, as each node is 'burned' it is connected to the 
GHZ state spread over the previously burned nodes. The 
process to do this is illustrated in Fig. [5] and consists of 




FIG. 10; After the singlets are formed we can repeatedly ex- 
tend a GHZ state from node A. This procedure uses the 
operation shown in Fig. [9] to add qubits to the GHZ state. 
The black squares depict qubits that are part of the GHZ 
state. Arrows represent a message to add nodes to the GHZ 
state along singlet paths. Each node keeps a record of the 
node from which it received this message from, symbolized 
here by a white dot. When a node cannot extend the GHZ 
state any further (highlighted by a dashed outline) it mea- 
sures its qubits in the X basis (open squares) and sends this 
information back towards A (thin arrows) along the route 
recorded. Certain nodes are selected beforehand not to per- 
form the measurement (here A and B) and these will form 
the resulting GHZ state. At each node the incoming data 
can be combined and sent back along one path if the routes 
branch. Once this data returns to ^ a phase operation can 
be performed on the qubit there to correct for any errors and 
the final GHZ state (here a singlet between A and B) will 
remain. 



joining the singlets partially contained in that node to 
the GHZ state. Within each node one qubit is left entan- 
gled with the GHZ state. After this operation has been 
run for a maximum of — 1 times all of the nodes in 
the cluster containing A have a qubit from a single GHZ 
state. 
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At each node a record is kept of the bond via which 
it has been included into the GHZ state. If there is a 
singlet between two nodes that are being burned then 
the singlet is ignored. Furthermore we add the rule that 
whenever a node can not extend the GHZ state anymore 
X basis measurements are performed along the recorded 
path back to A. This removes a qubit from the GHZ state 
but introduces a phase error in the remaining GHZ state 
depending on the outcomes of the measurement. The in- 
formation about these measurement outcomes has to be 
sent back along the path to A. Whenever the route back 
branches, the measurement outcome is sent in one way 
and a message corresponding to 'no phase error occur- 
rence' is sent to the others. At each node the returning 
process is paused until all of the bonds it sent a burning 
signal to provide it with the phase information. At nodes 
A and B we do not perform the X measurement. Finally 
after A receives all of the phase information a phase cor- 
rection can be performed and we obtain a singlet between 
A and B. In Fig. [in] an example is given to illustrate the 
protocol. 
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V. QUANTUM PREPROCESSING 

Despite being a very effective method, it is known that 
CEP in a network of pure states can be improved by cer- 
tain quantum 'pre-processing' strategies, and therefore 
CEP is not optimal [H, [11, M, [13, 111, [ll- In the fol- 
lowing we show that this is also the case in mixed-state 
networks. 



A. Swapping procedure 

To start with we generalize the swapping arrangement 
shown in Fig.[2]previously studied for pure states [13, [3 ■ 
In this arrangement we have two 2-qubit states that both 
have a qubit in a common node. If the two states are 
pure states \a) and |/3), with a > 1/2 and (3 > 1/2, we 
can obtain a singlet by swapping and then converting the 
resulting pure state into a singlet with a total probability 
of 2 min((l — a), (1— /3)) which turns out to be the optimal 
probability. Particularly CEP, which consists here of the 
Procrustean method followed by entanglement swapping, 
always has a smaller SCP of 4(1 — q;)(1 — Note that 
the optimal probability is equal to that of converting the 
least entangled of the two bonds into a singlet using the 
Procrustean scheme [l2j . 

To generalize this to mixed states we must consider 
double-edged bonds, each consisting of two PMSs, as il- 
lustrated in Fig. Ill[ since otherwise singlet generation 
would not be possible. Introducing more than one edge 
between the nodes allows us to concentrate the entangle- 
ment at different stages which gives rise to three different 
possibilities: 

I CEP - As previously described, the bonds arc converted 



FIG. 11: Three methods can be used to generate a singlet 
between two nodes A and C via an intermediate node B in 
an arrangement with two edges per bond. The thick black 
lines indicate pure but not maximally entangled states. 

to singlets and then swapping is performed over the 
resulting states. 

II Direct swapping - This applies entanglement swapping 

twice and then the resulting states are converted 
into a singlet. 

III Hybrid swapping - Here we distill a state of higher 

entanglement in each bond (but not necessarily 
a singlet) leading, if successful, to a single (par- 
tially) entangled pure state in each bond. This is 
followed by entanglement swapping and the Pro- 
crustean scheme to create a singlet. 

Each of these possibilities uses the swapping operation at 
different stages as illustrated in Fig. [TT] The exact im- 
plementations for the procedures depend on the types of 
states used. We will first apply each of them on a network 
of pure states and compare the SCPs. We then generalize 
to PMSs and show that direct and hybrid swapping can 
outperform CEP. 



1. Pure states 

If we start with bonds made of pure states |a) and |/3) 
we must have a way to convert each bond into a singlet 
in order to apply CEP (I). The method and highest possi- 
ble probability to accomplish this are given by Majoriza- 
tion [13 with a probability p = min(l, 2(l-a/3)) [ll,[ll]. 
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FIG. 12; Comparison of the three methods described in the 
text for creating a singlet between nodes A and C in the 
setup shown in Fig. [TT] for pure states. Shown are the success 
probabihties if the bonds are made up initially of the states 
I a) and |l/2) for CEP (solid line), direct swapping (dotted 
line) and hybrid swapping (dashed line). 



CEP applies this operation on each bond and if both 
bonds are converted into singlets swapping can be per- 
formed and the operation is a success. Therefore CEP 
succeeds with a probability (min(l, 2(1 — q;/3)))^. 

Our second method, direct swapping(II), is simply the 
application of the procedure for bonds containing one 
edge twice. If either generates a singlet the procedure 
succeeds. This gives a SCP of 1 - (1 - 2(1 - a))(l - 2(1 - 
/3)). There are adjustments we could make, for example 
use the results of Majorization to convert both of the 
states into a singlet with the highest possible probability, 
however all of these have a smaller SCP than CEP for a 
range of parameters. 

Finally, the hybrid swapping(III) method concentrates 
each bond to one pure state, | max(l/2, a/?)), with cer- 
tainty. This concentration procedure is also found using 
results from Majorization theory [s^]. Afterwards there 
is one pure state in each bond, as discussed previously, 
and we can then perform the strategy with optimal suc- 
cess probability min(l, 2(1 — a/3)), i.e. swapping over the 
pure states followed by the Procrustean method. We can 
actually consider the setup as a bipartite system between 
A and BC. The Majorization results then give the best 
possible probability of generating a maximally entangled 
2-qubit state between these systems as min(l, 2(1 — a/3)) 
which means that it must be the highest possible proba- 
bility for any method to succeed. 

Figure [T2l shows the probabilities in all three cases and 
we can sec that CEP is outperformed for a vast range 
of parameters by both other strategies. In hybrid swap- 
ping (III) , we have used multi-edged bonds to create pure 
states with the highest probability before applying en- 
tanglement swapping. We will refer to all strategies that 
have this property as 'hybrid'. This probability is unity 
for initial pure states but for mixed states the initial con- 
version of bonds to pure states is probabilistic, so when 
the conversion fails the bond is destroyed. 



S. Purifiable Mixed States 

We will now investigate if similar improvements can be 
obtained with PMSs, i.e. if the bonds between the nodes 
are composed of p(a. A) and p{P,h'). Again we will see 
that hybrid swapping provides the highest SCP. 

I CEP 

The classical percolation scheme involves perform- 
ing a PCM described in Sec. HIl followed by the pro- 
crustean protocol on both bonds and each succeeds 
with a probability given by Eq. ([5]) which simplifies 
to 

p,onv = 2Ai^min(a(l - /3),/3(l - a)). (25) 

To perform a swapping operation yielding a singlet, 
between nodes A and C we must succeed for both 
bonds which gives the total chance of success 

PC£P = (2Azymin(a(l-/3),/3(l-a)))2, (26) 

by simply squaring Eq. ((25)) . In this case the swap- 
ping operation is the final step of the protocol. 

II Direct swapping 

In our 2-edged setup we perform the swapping op- 
eration introduced in Sec. |TT] twice and there are 
two choices to do this if the states are not identical. 
Either we perform the swapping over the identical 
states p(a, X)® p{a, A) or we perform the operation 
on the states p(a. A) (8) p(/?, i^)- When we swap over 
identical states we obtain the state 

''( i-2:+2,^ '^'"-^''+^°'')- 

together with a further state where (3 is replacing a 
and ly is replacing A. Note that Eq. ((27| is obtained 
by setting 7 = 5 = and a = A = i/ in Eq. pT|) . 
This pair of states can then be transformed into a 
singlet with a probability 

Pd' = 2X^1^^ min(a2(i _ /3)2, /?2(i _ af) (28) 

which is calculated using Eq. In the case 

where we swap over non-identical states p{a, A) (g) 
p{(3, v) we obtain two states of the form (fTT|) with 
7 = (5 = 0. These can be converted into a singlet 
with probability 

Pd^2\^v'^al3{l-a){l- (3). (29) 

This is always larger than pd* and thus swapping 
with non-identical states should be preferred. 

III Hybrid swapping 

The hybrid method requires a concentration proce- 
dure to be performed (yielding a single pure state 
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FIG. 13: Success probability to generate a singlet between 
the end nodes of the swapping setup shown in Fig. [11] for the 
classical scheme (solid line), direct swapping (dotted line), 
and the hybrid scheme (dashed line). Each bond initially 
contains the states p(a, A) and p(l/2,A). We have indicated 
the percolation threshold of a face-centered cubic network. 



in each bond) which is given here by PCM. How- 
ever, in contrast to the pure state case discussed 
above, if a = /3 wc obtain singlets (in which case 
the method is identical to CEP) and, generally, the 
operation succeeds with a finite probability given 
by Eq. (HI). For non- identical PMSs PCM yields 
two non-maximally entangled pure states which are 
then used for entanglement swapping followed by 
the Procrustean method. The probability of suc- 
ceeding in converting both of the bonds to pure 
states is 

p2 = A2i/2(a(l-/?) + /3(l-a))2. (30) 
These pure states have largest Schmidt coefficient 
^ max(a(l-/3),^(l-a)) 



(a(l-/3)+/3(l-a)) 



(31) 



So, by using the SCP in single edged swapping with 
pure states we find that we can convert this pair of 
states into a singlet between the end nodes with 
probability 

^^ o min(a(l-/?),/3(l-a)) 

= ' (a(l-/3)+/3(l-a)) ' ^''^ 

Hence, the overall probability of succeeding with 
this scheme is 



Ph ^2\'^i^^[a{l - (3) + /3{1 - a)] 
X min[a(l — /?), /3(1 — a)]. 



(33) 



If we compare the success probability of direct swap- 
ping, pd, to the probability of success in the classical 
percolation scheme, pcep, it can be seen that classical 
percolation is more likely to succeed in producing a sin- 
glet if 

2min(a(l-/3),/3(l-a)) > max(a(l-/?), /?(l-a)). (34) 




FIG. 14: Illustration of entanglement percolation in a 3D 
network. The circles represent nodes containing qubits and 
the lines represent bonds containing pairs of two-qubit entan- 
gled states (the edges are not shown). The 3D network can 
be transformed into a Face-centered cubic network by per- 
forming the swapping operations (see Fig. [9]) over the smaller 
nodes. For some bond parameters the hybrid scheme allows 
percolation to occur where classical percolation fails. 



But the ratio of the success probability for the classical 
scheme against the hybrid protocol, ph, is 



PCEP ^ 2min(a(l-/3),/3(l~a)) 
Ph (ail - 13) + 13(1 - a)) 



(35) 



Whenever a f3 this is less than one and there is an 
improvement over the classical percolation scheme. Fur- 
thermore, the hybrid scheme is more likely to succeed 
than direct swapping. In Fig. [T3] we compare the proba- 
bilities of success for all schemes. As can be seen, hybrid 
swapping leads to the highest success probability. 

Hybrid swapping can be used in sections of larger net- 
works to allow percolation to take place. A simple exam- 
ple is a face-centered cubic (FCC) network, where every 
bond is split into two 2-cdged bonds (sec Fig. [T4|). When 
the above schemes are applied at the nodes linking two 
2-edged bonds the FCC network is recovered. Percola- 
tion is possible in these 3D networks with a threshold 
of approximately « 0.12. Since the classical scheme al- 
ways gives a smaller success probability than the hybrid 
scheme there are cases where the hybrid scheme allows 
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FIG. 15: Application of the hybrid scheme in a square net- 
work. This involves transforming the PMSs into pure states 
probabilistically and then applying a suitable pure state pro- 
cedure (see text). In the case shown all of the conversions are 
successful. When this happens a swapping operation can be 
performed and the resulting states distilled into a singlet. 

the percolation thresliold to be exceeded but the classical 
scheme does not (see Fig. [T^ . 

B. Square Protocol 
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FIG. 16: Comparison of singlet conversion probabilities for 
the different strategies in the square configuration, i.e. pcEP 
(solid line) and psq (dashed line) for X — v = 0.98, /3 = 0.5. 
We have also indicated the percolation threshold for a trian- 
gular network. 




CEP can also be improved on by using the hybrid strat- 
egy in a 2D square network, as shown in Fig. 1151 Each 
bond is converted into a pure state, |d) = %/d|00) + 
Vl " Q!|ll), by using PCM which is successful with a 
probability Pc on each bond. If this yields only two 
states \a) having a common node (B or C), entanglement 
swapping can be performed followed by the procrustean 
scheme. If all four PCMs succeed the resulting states can 
be connected (e.g. at nodes B and C) via a slightly mod- 
ified version of entanglement swapping, the so-called XZ- 
swapping . For this swapping operation the Bell mea- 
surement that usually has both qubits measured in the Z 
basis now measures one in the X basis. After this mea- 
surement unitaries are again applied to return the state 
into Schmidt form. The results of the Bell measurement 
have an equal probability, p„i = 1/4, for all outcomes m. 
Performing this operation twice on the square leads to 
two pure states (between A and D) of the form |a), with 
5 = (1 + y/1 - 16a2(l _ d)2)/2. These can be distilled 
into a singlet with probability min[l,2(l — a^)] by us- 
ing the protocol based on Majorization [30|. The overall 
chance of succeeding in generating a singlet is then given 
by 

Ps, = ApI{1-pI)[1 - a) +p^min(l,2(l - d^)). (36) 

When attempting to accomplish the same scheme using 
CEP we succeed with a probability ofp(7£;p = l"(l^ 
PcEp)^ which can be significantly smaller than Eq. (j36p . 
as shown in Fig. [TBI 



FIG. 17: By applying the square protocol on the shaded re- 
gions a triangular network of randomly distributed singlets 
(dashed lines in the right figure) is recovered. In the left fig- 
ure the nodes (circles) are linked by bonds (solid lines) each 
containing two edges (not shown). 



Again, this improved strategy may enable an infinite 
cluster to form when applied to larger networks. An 
example is shown in Fig. 1171 Here the square proto- 
col recovers a triangular lattice. If the conversion of the 
squares into singlets succeeds with a probability exceed- 
ing the percolation threshold an infinite cluster forms. 
In Fig. [16] it can be seen that the hybrid scheme exceeds 
the threshold for a triangular lattice in cases where CEP 
does not. 



C. Hierarchical Networks 

Small networks like the square configuration discussed 
above can be extended to larger networks in an iterative 
fashion. Networks formed in this way from pure states 
were considered in 



15[ . Again the probability of success- 



fully creating a singlet was shown to be larger when quan- 
tum strategies were used instead of CEP. However, the 
scheme with the highest probability is still unknown for 
these 'hierarchical' networks. Here we will consider two 
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different hierarchical networks with two edges per bond. 
Each of these contains the square network at some iter- 
ation level. We determine the SCP when using CEP in 
both cases which we then compare to the hybrid strategy. 
As it turns out, the hybrid scheme outperforms CEP. 

The first hierarchical network we consider is based on 
the 'Diamond' lattice, which at each stage replaces its 
bonds by the square network. The geometry for the first 
three iterations is shown in Fig. [THl The aim is to create 
a singlet between A and B and if we apply CEP the 
probability of succeeding at each level is given by the 
iterative formula 

pmarnond ^ ;l - (1 - p2_^^2^ (37^ 

starting with pf ^amond ^ ^^^^^^ 

The second hierarchical network we consider is the 
'Tree' network which is again built on the square con- 
figuration. For these networks an iteration is formed by 
creating two copies of the previous iteration and linking 
the bottom-left and top-right corner of the square to sep- 
arate nodes A and B as shown in Fig. [191 Again, we wish 
to generate a singlet between the opposite corner nodes 
{A and B) and CEP generates a singlet with a probability 

Pr'' = l-il-Pr-lplnvy, (38) 

where po = 1. 

Now we wish to see whether the hybrid scheme gives 
a larger SCP in these networks. Once again, the hy- 
brid scheme we consider starts by converting all of the 
bonds into identical non-maximally entangled pure states 
probabilistically. If the conversion fails on a bond then 
the bond is destroyed. This results in a network con- 
taining random pure state bonds. Each of these bonds 



FIG. 19: First three iterations of the tree lattice. Each it- 
eration is given by repeating the previous lattice twice and 
linking the pair of previous endpoints at new endpoints. We 
aim to create a singlet between nodes A and B. 



contains one edge. Ideally we would then apply a pure 
state protocol yielding the highest SCP between the in- 
tended nodes, however, this protocol is not known in the 
general case |15| . Instead we apply a procedure which 
performs XZ-swapping in cases when two bonds each 
have a qubit in the same node (except if these nodes 
are A or B). However, we also distill pure states into 
states with more entanglement whenever two edges form 
between two nodes and before performing further swap- 
ping. Finally, once one state is obtained between A and 
B, the Procrustean procedure is used to create a singlet. 

We applied this protocol to the hierarchical diamond 
and tree networks. For the second and third iterations 
of the diamond lattice the probabilities of creating a sin- 
glet are given in Fig. [201 together with the probabilities 
using CEP. This comparison was also made for the first, 
second and third iteration of the tree network and the re- 
sults are shown in Fig. [^l] These examples all illustrate 
an improvement in the probability of forming a singlet 
when using the hybrid method rather than classical per- 
colation. 



VI. CONCLUSION 

We have demonstrated that within lattice networks, 
where the nodes are connected by multiple bipartite 
mixed states, percolation strategies can be applied for 
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FIG. 20: Probability of succeeding in generating a singlet 
between the endpoints of a diamond lattice for the 2nd and 
3rd iterations (dashed lines). These give higher probabilities 
than the classical protocol (solid lines). The bonds contain 
two edges with parameters \ = v = 0.9, p — 0.5. 




FIG. 21: Probability of succeeding in generating a singlet 
between the endpoints of the 1st, 2nd and 3rd iterations of 
the tree lattice (dashed lines). These also outperform the 
classical protocol (solid lines). The bonds contain two edges 
with parameters \ — v = 0.9, /3 = 0.5. 



distributing entanglement. This is reliant on the states 
being PMSs, which are known from the DLCZ repeater 
scheme and arise as a result of amplitude damping. To 
show this we have introduced some new purification pro- 
tocols designed to maintain the form of these states or 
generate singlets. Like in the pure state case, a higher 
probability of distributing a singlet can be obtained, 
when the states in a bond are not identical. The ques- 
tion of whether quantum strategies can outperform CEP 
when each edge in a bond is identical is still open. Since 
we have shown that classical entanglement percolation is 
only possible for a specific class of bipartite states, en- 
tanglement distribution in a network which is subject to 
more general forms of noise needs to make use of other 
methods. The development of these methods is one of 
the most important goals for future work. These will 
not produce perfectly entangled states, however, the re- 
sulting state fidelity may be independent of distance and 
sufficient for purification. An example of such a strategy 



is given in Ref. [3l| for a bit-flip noise model. Progress in 
this direction has also been accomplished by gene rating 
3D thermal cluster states using Werner states [3^. Issj l . 
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APPENDIX A: PROOF OF SINGLET 
DISTILLATION REQUIREMENT 

In this appendix we give a concise proof for a neces- 
sary and sufficient condition to be able to create a singlet 
out of entangled mixed states using LOCC. We allow the 
states to be arbitrary bipartite states which are shared 
between two nodes and all operations are LOCC. A sim- 
ilar proof, but partly restricted to identical states, was 
given in |2l| . 

Lemma A.l. If a quantum state Pat can be distilled to 
a pure states, \^!) , then any state with the same range 
R{pab) is o,lso distillable to this state with non-zero prob- 
ability. 

Proof. The general form of the state is 

N 

Pab = ^Pt\tpt){'^t\, (Al) 
i=l 

with Pi > 0, X^jPj = 1 and e Ha® Hb. If the 
state is distillable to a pure state, j^*), there exist linear 

operators Ma and Nb, with MaM\ < I,NbN'^b ^ 
such that 

Ma ® NBPabM\ ®nI= pl^f) (*| (A2) 
^Ma®Nb\A){HmI®nIcx\^){^\ (A3) 

or 

MA®NB\^^)mM]^®Nl^O. (A4) 
This can be summarized as 

^ Ma®Nb\A) = q^\'^'), (A5) 

where at least one qi is non-zero as otherwise the oper- 
ator fails to distill pat- If this condition is satisfied the 
operation distills the mixed state into \'^). Now given 
another state p~ab with the same range as pab- We have 
that 

M 

P~ab^'^Pi\i'i){i^z\ (A6) 
i=l 
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with \'4)i) = Ylj=o0.t,j\'4>]) and = Y.f^ah.j\i'i) ■ This 
then gives 

AT 

Ma ® NB\i>i) = Ma ®NbY^ a,;^ |i/'j> (A7) 

JV 

= ^a.,,g,|*)=g~|vl/). (A8) 

The value of one must be non-zero as otherwise all qi 
are zero and this contradicts the fact that the operator 
distills pab- Hence the protocol also distills pab- □ 

Lemma A. 2. For n 2-qubit states to be distillable into 
a pure singlet at least two 2-qubit states cannot have a 
range spanned by product states. 

Proof. If a 2-qubit state has a range that can be spanned 
by product states then a separable state with this range 
exists. If there are n states each with a range spanned 
by product states the system state would have a range 
equivalent to a separable state formed by all of these 2- 
qubit separable states. Since it is impossible to distill a 
pure entangled state from any separable state it is impos- 
sible to distill a pure entangled state from n two qubit 
states each with a range spanned by product states. Sim- 
ilarly, if one of the 2-qubit states does not have a range 
spanned by product states, but all of the other states do, 
the range is equivalent to the range formed from a sep- 
arable state and one mixed 2-qubit state. This can not 
be distilled into a pure singlet as it would contradict the 
result in [22|. Hence at least two states can not have a 
range spanned by product states to be able to distill a 
pure entangled state. □ 

We now need to look at the two qubit states that sat- 
isfy this property. The states with rank one are already 
pure and if they have rank four the range can be spanned 
by product states. Similarly, if the state has rank three 
it can also be spanned by product states. This can be 
seen by considering the subspace orthogonal to a gen- 
eral state Y^|00) + y/1 — ajll). This space is spanned 
by {|01), |10), (^/T^|0) - V^|1))(|0) + |1]))/V2} and 
these are all product states. The last states to con- 
sider are those of rank 2, which fall into two cate- 
gories [sj]. The range is either spanned by product 
states {|00), (7A|0) -^/T ^|l))(V7Z|0 ) + VT^|1))} or 

{|00), (Va|01)-hv^|10))-FVl - - a|00))}. Hence only 
states that have a range containing one product state 
are the mixed states satisfying the condition. All mixed 
rank two states of two qubits can be considered to be the 
mixed state formed by tracing out a third qubit from a 
pure three qubit system. The classifications of these 3 
qubit systems is given in (ssl . [36l . [37j and for the range of 
the mixed system to contain one product state the three 
qubit state belongs to the W class. This class can always 
be written as \/A|$)|l) + Vl - A|00)|0) with 

1$) = V^IOI) + ^/p\10) + V7|00), a + + 7 = 1- (A9) 



By tracing out one qubit and using local operations the 
2-qubit state that can not be spanned by product states 
has the form 

p = A|V)(V'l + (l-A)|01)(01|), (AlO) 

where 

IV-) = V^|00) + V^|ll) + V7|01>, a + /3 + 7= 1. (All) 

So the only states that can be purified into a perfect 
singlet, given finite copies, are of this form. 

If there are two states of this form we know that the 
system is distillable since the procedure given in Sec. |TT] 
succeeds in the distillation. 



APPENDIX B: THE DISTILLABLE SUBSPACE 
SCHEME 

To extend the DSS scheme to n PMSs, p{a,X), we 
first need to describe the 2" non-zero eigenvalues and 
their eigenvectors. These correspond to different combi- 
nations of n — I \a) terms and / |01) terms. Then taking 
the decimal representation of the local states we can la- 
bel each of these eigenvectors by the decimal difference 
between the values at each location. This difference y in 
binary gives the location of the |01) terms. For example, 
in the case of two identical PMSs these are 



00 : 


\a)\a), 


10 : 


m\a), 


01 : 


|«)|01), 


11 : 


|01)|01) 



(Bl) 



and y takes all of the values from to 2" — 1. Now 
we define m{x) to be the number of Is in the binary 
representation of x and Ty = {x : x^y = 0, < x < 
2", X eN} ('^' is the bitwise AND operation). 

Then / = m(y) and all of the terms in a non-zero eigen- 
vector are of the form 

y^a"-™W-'(l - a)™(^)|x)l|a; + ?/)|, (B2) 

XGTy 

with eigenvalue A"~'(l — A)'. 

From this structure we can project out an entangled 
state if we measure the operator (|c)^(c|^ -I- |<i)54(d|l) at 
A and then (|c + y)%{c + y\% + \d + y)%{d + y\%) at B, 
when c & Ty, d & Ty, d > c and as long as there are 
no other terms of the form |c)^|(i + j/)^ -I- |(i)54|c + y)%, 
+ y)^ O'" + in non-zero eigenstate. 

The term |c)^|d -|- y)% + \d)'A\c -\- y)% can not appear 
in one eigenstate since all of the terms must have the 
same y value and this would require c to be equal to d. 
The state |c)j4|(i -I- y)g lies in one if and only if 3w e N, 
< t« < 2" — 1 such that c + w = d + y and c & Tw 
Similarly for |c?)j4|c -|- y)g but this case can not occur 
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since c G Ty and d G Ty means that w > y and d + w > 
c + y. If we assume that 3w gN, 0<w<2" — 1 such 
that c + w = d + y and c G T^j this would mean that 
d = c + k and that c^/c = for some A: > 0. Both of 
these resuhs then give that k^y = and w ~ k + y. So, 
ifw = y + k^d + y~c such that c^fc = can not be 
satisfied we create a maximahy entangled state. 

Now we have a choice of ways of creating these mea- 
surements. One particular way involves the definition of 
sets Sk = {x : m{x) = fc, < x < 2", a:: € N} and 
Ja,b = {x- a;^(a OR 6) = 0, < X < 2", a; € N}. Then 
the protocol consists of performing a POVM Pk.a.b = 



Cfe(|a>l 



l^)l(^ll) location A with a, 6 G Sk, 



a ^ b and < A; < n. For k = and n we define 
Pk^a,b = 12*^ — 1)1(2'^ — 1|^ and when these outcomes 
occur the procedure has failed. Here Cfc is a factor to 
ensure that 



k—n 

E 

a,beSk,k=0 



Pi 



k.a.b 



= I. 



(B3) 



With this outcome at location A another POVM is done 
at location B given by the operators = Ck{\a + 
d)%{a + d\% + \b + d)%{b + d\%) {d e Ja,b) and F = 
I — J2dGJ b ^<^' outcome here is F the protocol 

has failed, otherwise we have obtained a maximally en- 
tangled state. This protocol works since a, b £ Ty for all 
y G Ja,b but there is no w = fc -I- y such that k^a — 
and b + y = a + y + k, since if there were we would have 
b = a + k but m(a + k) ^ mib). 

The probability of succeeding is given by Eq. (fT6|) 
which comes from considering a particular eigenstate 
with parameter y. In this eigenstate there are 

^^^fn-rn(y)\ Un-rn{y)\\^^ ^^^^ 



different pairing terms, |a)l|a -|- y)^ -I- & + with 
a,b G Sk- The number of possible measured operators 
from this eigenstate is given by 



n — m{y) 
k 



- 1 



(B5) 



Note that the pairings in the eigenstate are twice as likely 
to occur than the ones with just an overlap and these have 
been counted twice in this sum. The probability that 
starting with an eigenstate (parameter y) we succeed is 
then 



2Ni 



1 



(B6) 



given that we have measured the operator Sk and the 
probability of this was 



Pk 



(B7) 



By summing over these we have, given we start with an 
eigenstate with m{y) = m, the probability of succeeding 
to be 



n — m — l 

E 



and these eigenstates occur with probability 



(B8) 



(B9) 



We have not counted = 0, n — m since they never con- 
tribute to the success probability. Then by summing over 
all of these we get the result in Eq. ((T| 
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